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Abstract 



We study the stabihty of the dislocation-free Bragg glass phase 
in a layered geometry consisting of coupled parallel planes of d=l+l 
vortex lines lying within each plane, in the presence of impurity dis- 
order. Using renormalization group, replica variational calculations 
and physical arguments we show that at temperatures T < Tq the 
3D Bragg glass phase is always stable for weak disorder. It under- 
goes a weakly first order transition into a decoupled 2D vortex glass 
upon increase of disorder. 
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An important problem is to understand how impurity disorder affects the 
translational long range order of periodic media such as the Abrikosov vortex 
lattice charge density waves (CDW) 0, Wigner crystals [0] and magnetic 
bubbles [Q. A purely elastic theory starting from the crystal, assuming the ab- 
sence of topological defects, such as dislocations, predicts an algebraic decay of 
translational order and that the resulting glass phase has divergent Bragg 
peaks [0 . This is in marked contrast with other models of glasses in superconduc- 
tors, such as the "vortex glass" or "hexatic glass" [050] where dislocations 



were argued to be present due to disorder beyond the Larkin length and trans- 
lational order decays rapidly. It was therefore claimed ^ that the results of the 
elastic theory are unstable to topological defects. However, in [H it was shown, 
using energy arguments, that as a consequence of the slow decay of correlations 
unbounded dislocations should not be present at equilibrium and at weak disor- 
der. The resulting phase, called "Bragg glass", is a glass with topological order. 
It has two important length scales, the Larkin- Ovchinnikov length Rc which char- 
acterizes pinning, and the translational order correlation length Ra- Ra can be 
much larger than Rc at low temperature if the scale rj a.t which disorder varies (i.e 
the superconducting coherence length ^) is small compared to the lattice spacing 
a. Since Rc is renormalized upwards by thermal fluctuations, one expects these 
two lengths to be of the same order of magnitude at higher temperatures and 
near melting. Such a phase is compatible with present decoration and neutron 
experiments [|ril-[T3[]. The random field 3D XY model being described by a very 
similar model, the same arguments in favor of a quasi-ordered vortex free phase 
apply at weak disorder. The existence of this topologically ordered phase in the 
3D XY model now finds some support from numerical simulations |T^. Also, 
detailed scaling arguments and a variational calculation supporting the existence 
of the Bragg glass were constructed recently in |T^ . 

In the domain of validity of the elastic theory in presence of disorder was 
estimated to be Ra ^ a. As we showed, this implies slow decay of correlations 
and, using energy arguments 0, the absence of dislocations. The Bragg glass is 
therefore a self-consistent solution until Ra ^ a. Beyond, the Bragg glass must 
undergo a transition into an amorphous state containing topological defects, upon 
increase of disorder or field. The nature of this state is unknown, various pos- 
sibilities are a strongly disordered vortex glass, an hexatic glass, or a strongly 
pinned liquid. A scenario suggested in [Q, is that this transition occurs around 
the field Htr of the observed tricritical point [T^ which separates, on the thermal 
melting line, first order melting at low field from second order vortex glass tran- 
sition at higher fields. From usual 3D melting considerations, one expects that as 
long as Rc ^ Ra ^ a the weakly disordered Bragg glass should melt through a 
first order transition upon raising temperature, while the amorphous state melts 
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through a second order transition. One can now find increasing experimental 
support for this scenario both in YBCO and BSCCO: the tricritical point can be 
lowered in field by controlled increase of disorder |T5[, the "second peak" around 



Htr appears as a sharp separation of two different types of behaviour in magnetic 
and transport properties [l^ll^l and neutrons |18[ also indicate a sudden loss of 
tranlational order at about Htr- 

The above self-consistency of the Bragg glass phase relies on arguments about 
dislocations. Also, to study analytically the transition between the two glass 
phases in c? = 3 is difficult since one must be able to describe the proliferation of 
topological defects. An interesting model geometry to study dislocation loops in 
presence of disorder was proposed by Kierfeld et al. [l^, who also derived more 
detailed energy arguments in favour of the existence of the Bragg glass up to 
Ra ~ a. The model consists in a model of coupled parallel layers of c? = 1 + 1 
lines constrained to lie within the planes and allows to describe some topological 
defects. It is relevant to the case of a magnetic field directed along the ab plane. 
This model was previously studied in the absence of disorder 



201, and found to 



exhibit a transition at Tc between a 3D coupled sohd for T < Tc and a decoupled 
high temperature 2D phase. Since the fiux lines are confined to the planes, an 
elastic description can still be used at weak disorder and therefore the methods 
of to describe the Bragg glass can be extended to this particular geometry. In 
the present paper we give a detailed solution of this model using both the rephca 
variational method and a renormalization group (RG) treatment. We show that 
at weak disorder the Bragg glass phase is stable, and that the criterion given in 
properly estimates the domain of stability of the model. Similar conclusions 
were reached by Kierfeld et al. using physical argument and a variational 
method. 

Interacting fiux lines confined to one plane in ci = 1 + 1, of average spacing a 
can be described in the continuum limit by introducing a smooth labelling phase 
field = 2nu/a, where u is the displacement field in the elastic limit. Since there 
are no topological defects within one plane the density of lines in each plane can 
be decomposed |]6|,|2l| as p{x) = po(l — Vn(a;) + X^p^o 6xp(ip(27rx + 0(a;)))) where 
p are integers. The model studied here consist of N planes of vortices coupled 
via their local density and in the presence of weak impurity disorder modelled by 
a random potential V{x) with correlations VkV-k = A^. The model is defined by 



the following random field Hamiltonian |22 



ij i 



X] 



flij COS( 
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with isotropic gaussian disorder (f{x)(j*{x') = iTgfjS^x — x') and 'Ua{x)rjp{x') = 
TAij6a/3S{x — x'). The original model in its simplest isotropic version, is defined 
by a smaller number of parameters: 

^ij^ = cSij ; fiij = ^{Si+ij + Si^ij) (2) 

4=^?%- ; Ai,=A6,, (3) 

= PqAp/T is proportional to the disorder in plane i with Fourier component 
close to 27ip/a, A = p^Ao/T is the long wavelength disorder coupling to slow 
variations of density [§. The extra terms present in (|l]), i.e the longer range 
couplings fiij, interplane disorder correlations gfj, long wavelength couplings and 
disorder correlators between planes, are generated by renormalization. They must 
thus be added to the general model which includes all most relevant terms allowed 
by symmetry. 

We start by applying the replica Gaussian Variational Method (GVM) to the 
replicated version of the starting Hamiltonian (^: 

H = jd'xY: \c{V4^? - /xcos(0^(x) - 

ai 

- E y V0^ - g co^mx) - 4>\{x)) (4) 

abi 

where a = 1, . . . , n is the replica index and one takes the limit n — at the end. 
For clarity have kept only the lowest harmonic g = g^, but the general case will be 
discussed below. Note that for a single harmonic model Rc ^ Ra 0- We study N 
coupled planes with ^ oo and use periodic boundary conditions. All planes are 
then equivalent and it is convenient to introduce the Fourier transform along z, 
rM) = 40'^(g,g.)exp[z(jg,)]. We denote !d\/{27:f by and / /„t/, rfg,/(27r) 
by where / is the interplane distance. 

One can approximate (Hf) by the variational Hamiltonian Hq = 



llq,q,Giq^qz)ab^''iQ^Qz)(p\-q,-qz)- The propagator G{q,q^)J is determined 
by optimizing the variational free energy F^ar = Fq + {H — i/o)//o- It is 
parametrized by a connected (thermal) part and a (disorder) self-energy part 
as = 5abG~^ — aab — Aq^ with '^ab = 0. The self-consistent saddle point 



equations [^] read: 



Ge{q, Qz) = l/icq" + 2p{l - cos(g,/))) (5) 



11 = p e 2^ii+i (ja^b = 2g e ^ " 



Bl\, = 2T f (1 - cos{qJ))Gaaiq, Qz) (6) 
B^^ = 2t[ Gaa{q,qz)-Ga,{q,q,) (7) 

•JCI.CIz 
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the latter being the intraplane and off-diagonal interplane mean squared phase 
ffuctuations, respectively. The quantity /t is the effective coupling between the 
planes and when it is non zero the problem is effectively 3D, whereas /i = is the 
signature of decoupling and corresponds to unbound dislocations being present 
between the planes. 

These equations have several types of solutions, corresponding to the three 
phases of the model: high-temperature (replica symmetric), 3D coupled elastic 
glass with full Replica Symmetry Breaking (RSB) and 2D decoupled glass (one 
step RSB). Each of these solutions is qualitatively similar to the ones obtained in 
0, for d = 3 and d = 2. Here however the model naturally exhibits transitions 
between these solutions. 

We first study the full RSB solution, as in 0. Parametrizing aab — (^(u) and 
j^ab _^ B(u), < M < 1, one has ct(-u) = 2ge~^^'^^"^^\ The inversion formula for 
hierarchical matrices gives: 

P r'^<= a'(v) 
BJu) = BJu,) + dv- ^ ^ 



q,qz -Ju 



icq^ + 2il{l-cosiqJ)) + [a]{u)y 



and Biiiuc) = 2T Jg ^^^cq"^ + 2/2(1 — cos(g2/)))^^ Differentiating the saddle point 
equation with respect to u, integrating over momenta and using [a]' = ua' one 
finds the full RSB solution for u < Uc- 

where the variables T = T/T^ and v = u/T have been introduced for convenience, 
Tc = Attc being the transition temperature of the 2D glass. For u > Uc one has 
[(t](m) = Si. It remains to find the breakpoint Uc, which is determined together 
with the effective coupling fl, by the equations a{uc) = 2ge~^"^^''^^'^ , and the 
equation for fl. Ba+i and the diagonal correlations Gaa{q,qz) can be computed 
using the inversion formula Gaa = Gc{l + Jq du [a]/v?{G~^ + [a]) and the above 
formulae. One finds two equations which determine Vc and fx: 

where A = cg^^^ ~ c(27r/a)^. The following functions have been defined, 
g{v) = (1 + t;)(^+ i)/^(l - i;)(i-^)/2t;-i and f{v) = g{v)(l){v) with 0(?;) = 
2we(^+^V(i-'')/(i+'^)-7(l + ^l-v"^). lYiom. the second equation in (0) one re- 
covers the decoupling transition at T^. for the pure system, and the corresponding 
value of fl for T < Tc using that f c — > when (7 — > and /(O) = e^. 
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For A = the equation for Vc is remarkably simple and reads gT/fi = 0(fc). 
In that case, and in the limit of very large cutoff, the transition from the 3D 
to the 2D vortex glass RSB solutions is continuous and happens when fc=l, i.e 
Uc = T. The transition line is thus determined by the equation: 

eT , , 

^^ = ^9 (11) 

Indeed when f c — >■ 1 one recovers smoothly the one step solution for the decoupled 
2D glass obtained in Ref. with Uc = T/Tc. 

Several effects (such as a finite cutoff, a small non zero A) transform this 
second order transition in a (weakly) first order transition (at weak disorder). 
When A > the above equation becomes: 



,5 A gT 
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<P{v,) {g{v,)r = - [^—^ (12) 

with 5 = AT/(Tc-T). 

The function in the left hand side reaches its maximum before f c = 1 and 
thus upon increasing T, the RSB solution ceases to exist before the point Vc = I 
is reached. Thus there is now a first order jump to the 2D decoupled phase. The 
value of Vc at the jump is given by f * = {1 — 6^)/{l + 6^). The equation of the 
transition line is then: 

The value of fi at the transition, i.e the size of the jump, is given by: 

(l-f)(l-5) .rl-f 
iL] - ^ ^S(l+T)-l /..N 

A J A{l + 5y ^ ^ 

The above result (|lTD for the transition from the 2D to the 3D glass can be 
recovered in a simple way by considering the characteristic lengths of the problem. 
In the 3D phase (at relatively large /i) one can expand the interplane coupling 
terms and obtain a manifold with bare elastic coefficients: cn = c(27r/a)^/Z 
isotropically in plane and Cz = fil{27i/a)'^ along z. The 3D translational order 
correlation length Ra along z was computed in as R^'^^ = a'^cnCz/n^A with 
here A = 2Tg/l. One then finds that the above transition corresponds to the 
ratio R^z/^ being a fixed number R^z/^ ^ e = 2.7... This can be viewed as a 
"Lindemann criterion" and is in reasonable agreement with the one obtained in 
1^ by a different variational method. The physical interpretation is thus that the 
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transition happens when the 3D translational order correlation length becomes of 
order the interplane spacing. Note that in general it should be the translational 
order length Ra which controls the transition and should enter in a Lindemann 
type criterion, and not the Larkin length R^. Although for a single cosine model, 
valid at higher temperature or ^ ~ a, one has Ra ~ Rc, when one treats model 
(1) keeping all harmonics in order to describe the low temperature region, one 
finds indeed that it is Ra which is the relevant length. 

To interpret further the above results in a simple way one can notice than in 
the equations for fi, relevant Qz momenta are the one close to the zone boundary 
[Qz ~ ±7r). The correlation function appearing in the equations for fl are therefore 
analogous to two dimensional ones, but for the presence of a mass term of order 
jl in the denominator. One can therefore crudely replace the equation for fl by 

where -B2d('") is the correlation in the purely two dimensional (disordered) system. 
On this expression it becomes clear that if correlation functions in 2D decay too 
rapidly as a function of distance there will be no solution of (|15]) for small fl, and 
the transition will be first order. To get a continuous transition one need that 
the right hand side of (|15D goes to zero with ft, at worst with a finite slope. This 
imposes at worst an algebraic decay of the correlation functions with an exponent 
smaller than one. This is the case for the pure system below T^, and for the 
disordered system with infinite cutoff and no A, where the exponent is frozen at 
the Tc value (i.e. one). In the presence of A > 0, or if one takes solutions given by 
the Cardy-Ostlund RG, for which correlation functions decay as e"^"'-'^^ , one gets 
a discontinuous transition, as confirmed below. A similar discontinuous transition 



was found from the effect of a finite cutoff in [|19[ . Based on (|T5|) one can construct 
an argument from the 2D translational order correlation length for a single plane. 
It was estimated in g as i?^^ = a(27rA/T5()i/2{i-T/Tc)_ rp^.^ j^^g^j^ should be 

compared with another 2D length for the pure system, i.e the scale below which 



the coupled system is still 2D, R2d = 2vryc//2 = a(/ie"^/A)^/^*^^~^/'^'=\ Remarkably 
the transition again occurs when these two lengths become comparable. 

We now turn to the Renormalization Group (RG) method, which can be 
applied to this model near the decoupling transition of the pure system, i.e near 
T = Tc = 4:71, for small disorder g and couplings between planes. We study 
the model (|lD. Near Tc one needs to keep only the lowest harmonic p = 1, 
higher harmonics being irrelevant. We have obtained the (replica symmetric) 
RG recursion relations for Kij, Aij, gij and fiij using Coulomb gas techniques 
and fermion methods. This is a generalization of the Cardy Ostlund recursion 
relations |^ to a set of coupled planes. For simplicity we give here only the 
results for two planes [^]. The results for planes will be presented in but 
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FIG. 1. Schematic phase diagram for two coupled planes in presence of disorder 



our preliminary studies show that they are qualitatively similar (see below). One 
defines, for two planes: K~j^ = (T/Tc)5ij + kij, kc = ku + ku, k = ku — ki2, 
= fJ'12 + 912, 9 = 9ii and 6 = An — A12. One finds the following RG flow 
equations upon a change of cutoff a — a e' 



; ^.^(J^-,?.) (16) 

—g = (^k + kc)9 + ^^9i2-9u~9'^ (17) 

^5-12 = (-25 + k + k^)gi2 + H9 - ^9129 (18) 

j^fi = 2{-6 + k)fi-g,29 (19) 

and j^kc = 0, J^An = 9'^/ 4. These equations are valid to second order in all 
coupling constants and in deviations from T^. The parameter kc = {T^ — T)/Tc 
is the true reduced temperature and is not renormalized. k controls the stiffness 
associated to the phase difference between the two planes, and is renormalized, 
indicating whether the planes are coupled or not. 

The previously known limiting cases are the following. In the absence of 
disorder, i.e 9 = 912 = 5 = 0, one simply recovers a sine-Gordon model for 
the relative phase 0i — 02^ there is a usual Kosterlitz-Thouless transition, with 
a separatrix k = —fi/2 between a coupled phase /i(/) 00 and a decoupled 
high-temperature phase (/u(/) —^0). In the absence of coupling between planes, 
i.e /i = 0,9i2 = 0, A12 = 0, A; = 0, there is a a 2D glass phase for kc > - 
characterized by the Cardy-Ostlund attractive fixed point 9* = kc and 5 (/)—*> 00 
- as well as a high temperature phase for kc < 0, where disorder is irrelevant and 
9{l) and 6{l) 6*. 

In presence of both disorder and interplane coupling, we find, via numerical 
solution of the equations that there are three very different regimes of the flow 
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with marked separatrices between them (we start from the initial condition k = 
kc,gi2 = 0, Ajj = corresponding to model (4)). There is a decoupled 2D glass 
regime (/i(/) 0,gi2{l) ^ 0,6 ^ oo, k{l) — > k*), a 3D coupled glass phase (all 
constants become large) and a high temperature phase gi2{l), g{l) 0, 

6{l) — > 6*,k{l) k*). These regimes are very likely to correspond to three 
different phases, the uncertainty on the exact position of the separatrices being 
quite small [^]. We thus arrive at the phase diagram of Fig. 1. In particular the 
phase transition between the 2D and 3D glass phases occurs approximately for: 

fi^C g (20) 

with C ~ 0.5 — 0.6, and only a weak dependence on temperature. This result 
is in good agreement with the result for the transition line from the GVM. One 
also notices that, when approaching the 2D-3D glass separatrix from the 2D 
glass phase, the interplane coupling always becomes vanishingly small at a fixed 
characteristic length l*{g), which diverges as g ^ but remains fixed close to the 
transition. We find numerically that l*{g) ~ l/{2kc) log(l/5') for small g and thus 
we identify this length with the 2D Larkin length. This may be interpreted as 
the 2D to 3D glass transition being first-order, also in agreement with the GVM. 
Finally, the study for N planes gives similar results for the phases. The separatrix 
remains roughly of the form (pO|), with C ~ 0.7 — 1.0, but deviations from the 



linear behaviour is observed p3[ at small g and near T^. 

In conclusion we have shown that in a layered geometry, the Bragg glass phase 
is stable at weak disorder. Upon increase of disorder it undergoes a transition into 
a topologically disordered glass. In this geometry this glass is in fact a decoupled 
2D glass. The transition is found to occur roughly when the translational order 
length Ra is of the order of the lattice spacing a as suggested in This model 
has some peculiarities that make generalization to real vortex lattices difficult. 
However, the fact that it confirms the existence of a stable Bragg glass is an 
additional indication, besides the energy arguments of 0, that the Bragg glass 
occurs in real vortex lattices. It also confirms that the criterion Ra ^ a may 
indeed be used as stability criterion for the Bragg glass. Our conclusions agree 
with the one of |]19| . 

We are grateful to T. Nattermann for interesting discussions which motivated 
this work. 
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